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Classical defects (monopoles, vortices, etc.) are a characteristic consequence of many phase 
transitions of quantum fields. Most likely these include transitions in the early universe and such 
defects would be expected to be present in the universe today. We continue our analysis of the onset 
of classical behaviour after a second-order phase transition in quantum field theory and show how 
£Nj defects appear after such transitions. 
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I. INTRODUCTION 



Because phase transitions take place in a finite time causality guarantees that, even for continuous transitions, 
correlation lengths remain finite. Order parameter fields become frustrated, and topological defects arise so as to 
■ reconcile field phases between different correlated regions . 

A huge variety of defects is possible, according to the complexity of the initial symmetry and its breaking [0. 
Monopoles can easily overwhelm the energy density of the universe, while cosmic strings (vortices) may be the source 
of high energy cosmic rays and lensing, as well as contributing to structure formation (although they are not now 
, thought to be its determining factor). More complex strings can exacerbate baryogenesis. 

As solutions to the classical field equations, these defects have non-perturbatively large energies, commensurate with 
the temperature scale at which the transition takes place. Thus, for example, cosmic strings produced in the early 
universe will be expected to have an energy/unit length (tension) a ~ (10 16 GeV) 2 , where 10 16 GeV is the estimated 
GUT scale. For this reason alone, defects are manifestly classical entities after the transition, whose evolution can be 
handled successfully through the solution of the classical equations that they satisfy (e.g. the Nambu-Goto action for 
tIh ' cosmic strings). Nonetheless, their origin in the beginnings of the phase transition that spawned them is quantum 
mechanical. There are further complications according as the symmetries are global or gauged Q, but the simplest 
pi 1 of all are topological global vortices (or strings) and it is these that we shall consider here. 

Having set up the model, the subsequent sections of this work concern the transition from quantum field theory to 
i-C ■ classical field theory. We know the mechanism for this, the decohering effects of the quantum environment with which 
the long- wavelength modes of the order parameter fields interact ||. Field or phase ordering and classical defect 
formation are controlled by long-wavelength modes. We consider these modes to constitute the open 'system' which 
undergoes quantum decoherence due to an 'environment' which consists of everything else (short-wavelength modes 
and all other fields interacting with the order parameter fields) . This permits the calculation of a decoherence time 
to after the onset of the transition, after which the long- wavelength modes behave classically, subject to generalised 
Langevin equations. At the same time, the Wigner functional can now be interpreted as a Fokker-Planck probability 
density and we can calculate correlation functions Given that the correlation functions are now defined in terms 
of classical probabilities, classical defects can be identified. 

The major characteristic of defects is their topological charge. We conclude by showing how localised topological 
charge precedes the appearance of defects, and gives an estimate of defect densities when they do appear. 

For an instantaneous quench the onset of classical behaviour has already been discussed by us in an earlier paper 
in these proceedings, henceforth known as I (see Q), and elsewhere ||, and we shall not repeat all the details here. 
What we shall do is elucidate those properties of the classicalisation due to the environment that are relevant to defect 
formation. 

We conclude by briefly comparing our approach to others for the production of defects. This work has yet to be 
completed, and will be continued elsewhere. 
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II. THE MODEL: COUPLING ORDER-PARAMETER FIELDS TO THE ENVIRONMENT 



The most immediate difference between the modei here and that of paper I |7|] lies in the replacement of the 
single real scalar order parameter field there by a complex field This is the simplest field permitting acceptable 
defects, in the form of global vortices. [The defects of a single scalar field with broken Zi symmetry, domain walls, 
would dominate the energy density of the universe if they were present.] As in I we take the $ field to interact 
with a collection of real scalar environmental fields Xn (n=l,...,N). It is no loss of generality to take all couplings to 
X-fields identical, and all Xn masses identical, to give a model with t/(l) x O(N) or 0(2) x O(N) symmetry broken to 
O(N). However, the O(N) symmetry of the model, while tactically useful, is accidental and purely for calculational 
convenience. We adopt a Cartesian field basis <& = {4>\ + i^>2)j\f2 because, in the linear regime at the early times 
that are most relevant to us, §\ and 02 behave independently. The action is 

S[$, X ] = SsystI*] + SenvM + $nt[$,X], (1) 

where we have made a provisional separation into system field and environment with which it interacts, where 
is the action for the <fi a fields (a = 1,2) (we assume summation over repeated index a), with fj? > 0, and 

ScnvM 

describe the environmental fields Xn and its interaction with them. These are taken to be weakly coupled, with 
1 l/VN ^ g ~ A. Meanwhile, for simplicity the Xn masses are fixed at m ~ fj,. 

The model has a continuous transition for the breaking of the 0(2) symmetry at the critical temperature T c where, 
in units in which ks — 1, T 2 = 0(fi 2 /y/NX) ^> fj 2 . We shall assume that the initial states of the system and 
environment are both thermal, at temperatures T = 0(T C ) > T c . In this way the <I>-field is peaked strongly in field 
space about the unstable maximum. On incorporating the hard thermal loop 'tadpole' diagrams of the x (and $) 
fields in the <& mass term leads to the effective action for $ quasiparticles, 

s^ st m = J d^x {1 wv a - \n\{T^\ - \m~\ 

where m|(To) = — fj, 2 (l — Tq /T 2 ) — M 2 > 0. As in I, we can take an initial factorised density matrix at temperature 
To of the form p[Tq] — p$[Xb]/3 x [To], where p$[To] is determined by the quadratic part of S , ^ f st [$] and /3 x [To] by 
<Senv[x]- That is, the many \n fields have a large effect on but the $-field has negligible effect on the Xn- This is 
crucial for the validity of our approximation. More details are given in I. 

In fact, with defects in mind, we need to perform a further factorisation to isolate the long-wavelength modes of the 
system field. As we stressed, the decohering agents are not only the environment fields Xn that have been introduced 
already but also the short-wavelength modes || of the self-interacting $-ficld. 

As in I, we split the fields initially as $ = <!>< + <!>> (and 4> a = 4> <a + 0> a ) where the system-field $< contains the 
modes with wavelengths longer than the critical value and the bath or environment- field <!>> contains wavelengths 
shorter than This is a natural division, since modes with wavelengths shorter than are stable, modes with 
longer wavelengths unstable. Further, the short-wavelengths explore the interior of classical vortices, irrelevant for 
counting them. However, once the power in the fluctuations has moved to long-wavelength modes, with wavenumber 
ko < fj,, the exact delineation is unimportant, and shorter wavelength modes can be included in the 'system'. With 
this in mind a more realistic decomposition into 'system' and 'environment' is not (]]]), but 

S[$, x] = 5 syst [$<] + 5 cnv [x, $>] + S lnt [$<,x, $>], 
with [i the demarcating momentum, where 

s syst [$<] = J d 4 x |ia^ <a 9^ <a + i M V< a - ^< a ) 2 } , 
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and 



^ 2 >a-^ 2 >a) 2 



Se„v[<I>>,x] = J ^{~^> a ^ >a -~ 

n=lV ^ J 

iV 

„=1 8 viv J 

- 2 / ^ ^<a( X ) < l > >b( X ) + 2(j) <a (x)(f> <h (x)(l) >a {x)(j) >h {x)}. 

All terms omitted in the expansion are not relevant for the one-loop calculations for the long-wavelength modes that 
we shall now consider. 

III. THE SIMPLEST TRANSITION: INSTANTANEOUS QUENCH 

For simplicity we repeat our assumptions in I, in adopting an instantaneous temperature quench from To to Tf = 
at time t = 0, in which m^(T) changes sign and magnitude instantly, concluding with the value m 2 ^ = —p 2 , t > 
(and beginning with the value m^(T ) = m 2 = 0(p 2 ), t < 0). We stress that is the renormalised mass, containing 
the temperature dependent interactions with all fields. 

A. The influence of the environment 

As we observed, the most visible signal that the transition has occurred will be the presence of topological defects, 
solutions to the classical equations of motion SS{&, x] — 0, with Xn = 0. These are global vortices in the field, around 
which the field phase 9 (3>(x) = h(x)e l8 ^) changes by 2n. Considered as tubes of 'false' vacuum, with cold thickness 
0(/Li _1 ), they have energy per unit length a — 0(p 2 /X) (up to multiplicative logarithmic terms ln(/i£dcf), where £def 
is vortex separation). In particular, the field $ vanishes along the vortex core. We can therefore use line zeroes to 
track classical defects Since classical defects are specific field profiles we need to be able to distinguish between 

different classical system-field configurations evolving after the transition. As a result, we are only interested in the 
field-configuration basis for the density matrix. 

Since the effects of (bosonic) environments are cumulative each contribution to the environment increases the 
diffusion term and thereby speeds up the onset of classical behaviour. Thus any part of the environment sets an upper 
bound on tjj. For large (but not infinite) N a one-loop approximation is sufficient for calculating tp due to x~ fields 
alone. It is convenient not to have complex arguments and we use <$> to denote the real Cartesian doublet (4>i, 4*2) and 
X to denote (xi,X2, ")Xn)> etc. At time t > the reduced density matrix p r [</><, , t] — (</><|p r (i)|</><) is now 

Pi#<i0<j*] = \ V X J V( t ) > /#<'0>>X> i4><,4>>,xA, (2) 

where p[<j><, </>>, XI <f><> 0>: X! t] is the full density matrix, £></>> = 2?</>>i2?</>>2 and T>x = rii ^Xn- 

The environment will have had the effect of making the system effectively classical once p T (t) is essentially diag- 
onal. This is very different from the late-time dephasing effects found in Refs. pO| , pH , which rely on time-averaged 
diagonalisation. More details are given in paper I and Jm]. Quantum interference can then be ignored and we obtain 
a classical probability distribution from the diagonal part of p T (t), or equivalently, by means of the reduced Wigner 
functional. For weak coupling (see I and 0) there will be no 'recoherence' at later times in which the sense of classical 
probability will be lost. 

The temporal evolution of p x [0< , 0< , t] is 



p T [cj>+,<p-,t} = / Xty+ / P0r < J I .[0+,0 f< ^|0+, ( / )i< ,t o ] p^+^io], 
where = T>cj>t 1 T>(f)t 2 , etc. and J r is the reduced evolution operator. 
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In order to estimate the functional integration which defines the reduced propagator, we perform a saddle point 
approximation 

JM^ATkA^^M « ex pi A[0+ cl ,0- cl ], 

where 0< cl is the solution of the equation of motion < \ < p+ = ^ ) - = with boundary conditions </>< cl (io) = </>j< an d 
^< cl (0 = It i s very difficult to solve this equation analytically. In order to make it tractable we assume that the 
system-field contains only one Fourier mode with k — fco for the reason indicated earlier, that the long-wavelength 
modes, for which | fco 1 2 < A* 2 ; increasingly bunch about a wave- number fco < A 1 which diminishes with time. However, 
unlike in I, we are interested in more than the fco = mode. 
For such small fco the classical solution is of the form 

= f a (s,t)cos(k .x), 

where / a (s,i) satisfies the boundary conditions / a (0,i) = </>i< a , and f a (t,t) — </>f< a . Qualitatively, / a (s,i) grows 
exponentially with s for t < t sp , and oscillates for t sp < s < t when t > t sp . For t < t sp we approximate it by 

sinh[<j (f-g)] sinh(o;os) 
7a(s,t) = 0i< a —, — - — + f<a ^-- 7 — — (3) 

where ujq — — fcp. 

We saw in I (see also that the linear approximation is reasonable until almost t ap , where the spinodal time t sp 
is defined as the time for which (| < I , <| 2 )t ~ i] 2 . That is, t sp is the time it takes for the field to populate the ground 
states of the model. As a result, t sp is given by 

exp[2/ii sp ] « O(^r). (4) 

The exponential factor in Eq.(^), as always, arises from the growth of the unstable long- wavelength modes. The factor 
X^T 1 comes from the coth(/3a>/2) factor that encodes the initial Boltzmann distribution at temperature To > T c . Thus, 

A^s P ~ ln (^==)- (5) 

As in paper I, it is sufficient to calculate the correction to the usual unitary evolution coming from the noise kernel. 
For clarity we drop the suffix f on the final state fields. If A = (| < I >+ | 2 — |$~| 2 )/2 for the final field configurations, 
then the master equation for p r ((f><, t) is 

^p r = (0+|[#,p r p<) - iVA 2 D(k a ,t) Pl . + ... (6) 

The time dependent diffusion coefficient D x (ko,t) that determines the effect of the \- fields on the onset of classical 
behaviour acquires a contribution D a (ko,t) from each field Xn, 



D n {k ,t) 



16N 



[ ds u(s,t) [ReG 2 ++ {2k :t- s) + 2RcG 2 ++ {Q;t- s)] . (7) 



where, for the case in hand of an instantaneous quench, u(s,t) — cosh 2 Lu (t — s) when t < t sp , and is an oscillatory 
function of time when t > t sp . The G++ are the long- wavelength correlation functions of the x-fields for the appropriate 
contours in the closed time-path. 

The contribution from the explicitly environmental ^-fields, D x (ko,t) = ^2 n D n (ko,t), takes the form 

a 2 T 2 

D x (ko,t) ~ ^^o exp[2wot], (8) 
p 6 

largely from the end-point behaviour at s = of the integral (0). For t > t sp the diffusion coefficient stops growing, 
and oscillates around D(ko,t = t sp ). 

The environment is also constituted by the short-wavelength modes $> of the self-interacting field. This gives an 
additional one-loop contribution to D(ko,t) with the same u(s) but a constructed from the short-wavelength 

modes of the $-field as it evolves from the top of the potential hill. Omitting this mode in its calculation just means 
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that any decoherence time to obtained from the x„ alone will be an upper bound on the true decoherence time. 
An estimate of the $-field contribution, based on the inclusion of tadpole diagrams alone, suggests that, with no 
l/N damping, the <!>> modes have the same effect on the dissipation, qualitatively, as all the environmental fields Xn 
put together. At an order of magnitude level we can ignore the <&> modes in the calculation of to, for which the 
%-fields alone give a strong bound. We cannot ignore them in the calculation of the defect density, since it is the 
coarse-graining of the order-parameter field that renders the line-zero density finite. 

We estimate to by solving (in the one- loop approximation) for the off-diagonal elements of the reduced density 
matrix 



4><A ~ P?[</>< ><£<;*] exp 



VT / ds D x (k ,s) 



(9) 



where is the solution of the unitary part of the master equation (i.e. without environment). In @ T — 0(p 4: ((j)S) z 



1^1 



<&^|)/2/z and S = (|$<| — \$> < \)/2p. V is understood as the minimal volume inside which there is no 



possibility for coherent superpositions of macroscopically distinguishable states for the field. We take this as 0(p~ 3 ) 
since p~ x is the thickness of an individual vortex. Inside this volume we do not discriminate between field amplitudes 
which differ by O(p), and therefore, as in I, take S 2 ~ 0(1). Similarly, we set cjr ~ 0(a/X), where A < a < 1 is to be 
determined self-consistently. 

The decoherence of the long- wavelength fco-mode by the environment occurs when the non-diagonal elements of the 
reduced density matrix are much smaller than the diagonal ones. In (^) this corresponds to when 



1 > VT / ds D x (k ,s) 



(10) 



since p"[<^<, 4 > < ; t] is increasingly independent of 5. See I for more details. Because the diagonalisation of p r {t) occurs 
in time as an exponential of an exponential, decoherence occurs extremely quickly, at time pv 



uj t D > ln(— -=). 

±c\/a 



(11) 



For u> ~ p, the value of a is determined as a ~ p/T c from the condition that, at time to, (\4'\ 2 )t ~ cvn 2 . Since 
a <C 1, in principle, the field has not diffused far from the top of the hill before it is behaving classically. 

It follows from Eq.([ll|) that the upper bound on to and, we assume, to itself, increases as k — > /i, although we 
need a better approximation to sec how, in detail. However, we stress that, as far as counting vortices is concerned, 
all that matters is how the power in the field fluctuations is distributed. The distance between defects is the relevant 
wavelength, and not defect size. With coq « \jl for the relevant k 2 = O^/to) it follows that 1 < pto < pt sp , with 



1 T c 
pt sp - pt D ~ - ln( — ) 
4 p 



(12) 



Whereas the environment is very effective at decohering adjacent field configurations, it has much less impact on 
the diagonal matrix elements 

Pt[tf><] = (0<|fr(i)l^<>, 



< at time t, once the theory is classical. This is 
<, and thereby p^,H to calculate the density of 



which give the relative probability that the field takes the value 
all that is needed to calculate equal-time correlation functions of 
classical defects by itemising their zero-field cores. 

We know that, for very early times, the Gaussian approximation for p r [4> <1 cf> < , t] is valid, although we should not 
interpret it as a probability then. We shall now argue that can we use the approximation Pt[4><] — p" [4>< , 4>< > t] j 
compatible with (|E)|) at least until time to- That is, although the environment is crucial in making p t off-diagonal, it 
has much less effect on the diagonal matrix elements, which are the ones we use for calculation. 



IV. CLASSICAL EQUATIONS 

Our analysis of the onset of decoherence shows that it makes little sense to talk about line zeroes of the field before 
to of Eq. (pi]) , since they would be expected to suffer from quantum interference, as well as having a density that is 
strongly dependent on the scale separating </>< from (j> > . 
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By the time t sp the Gaussian approximation has broken down, and the Goldstone field phase 9 will have decoupled 
from the heavy h-mode. From here onwards it is the massless Goldstone modes whose causal propagation controls 
field ordering. This is the last ingredient in turning a line-zero into a classical vortex. 

This requires both classical probabilities and classical equations. Let us consider probabilities first. The reduced 
Wigner functional is defined as 



V V< e"<"< (^<-ry < |p r (t)|0 <+??< > 



whereby 



P[</><] t = / VlT < W t [</><, 7T<] = (0<|p r (i)|<M 



(13) 



(14) 



is the probability density for field configurations. 

For t > tjj of (pi]), W t [4><i 7r <] is positive, at least for long-wavelengths, and can be identified with the Fokker- 
Planck probability distribution function P t FP [<j) < , 7r<] , from which P[<^<] t can be equally identified as the Fokker-Planck 
probability 



P 



FPr 



} <\t 



J 2^<P t FP [<P<, 7T<]. 



(15) 



Suppose we can calculate P[4><]t for times t >to- This permits us to calculate the equal time n-point correlation 
functions 



G^ fc .. c ( Xl ,..,x n ,t) = J D0 < if p [0 < ]^ <a (x 1 )..^ <c (x ll ). 



(16) 



As we know [g_3|jlj], equal-time correlators are all we need to calculate densities of zeroes (line-zeroes, etc.). Moreover, 
the dominance of P[(f><\t by long- wavelength modes permits the adoption of a single decoherence time tjj from Eq. ([l2|) . 
However, if all we are going to use is P[0<]t, the diagonal matrix element of pt, there is no real need to construct the 
Wigner functional. We can just do a calculation of P[0<]t from the start, along the lines of Ref. 15|. 

To see the appearance of individual defects is more difficult, and requires classical stochastic (Langevin) equations. 
We have yet to complete the analysis. However, as a starting point let us pretend that the only environment is the x 
fields. In the usual fashion, one can regard the imaginary part of 5 A as coming from a single noise source £,(x), with 
a Gaussian functional probability distribution given by |p|, [l6|jl7| | 



T>[£\ = iV 5 exp 



g 2 N 



(17) 



where iVg is a normalization factor, and N(x — y) cx ReG^_ + (x — y) for the single x-loop. Indeed, we can write the 
imaginary part of the influence action as a functional integral over the Gaussian field £,(x), 



2?£P[£] exp 



d A x- i\ A(x)£(x) 



exp 



i Id x Id y 



A(x) g 2 N(x iy ) A(y) 



Therefore, the coarse-grained effective action (see I for definitions) can be rewritten as 

A[<t> + ,4>-\ = J VtP[£] exp | zSeff [</>+, 



where 



Taking the functional variation 



d A x 



A(x)£(x) 



(18) 



(19) 



(20) 



G 



gives the "semiclassical Langevin equation" for the system-field ^j,^6|JT^] (up to factors 0(1)) 

□ a (x) - /i 2 a + \fa(x)<fc(x) + g 2 ^{x) [ d A y K{x, y) (&{y) = <f> a (^{x), (21) 



where we have assumed only quadratic interactions, and K is the (assumed) common mass retarded loop, arising 
from the real part of 5 A. jl and A are the 'renormalised' constants by virtue of the coupling with the environment. 
To estimate the size of g 2 cf> a (x) J d 4 y K(x, y) <p^(y) in comparison with the /z 2 </> a term, we note first that 

K oc TV - T C V - ///ViVA (22) 
for high T. Further, with relevant distance scales no larger than /i -1 then, crudely, at time t sp , 

\g 2 J d 4 yK(x lV ) cf>l(y)\ 
<g 2 f d A y \K(x,y)\ cb 2 (t sp ) 



-A 2 ^t sp ^(t sp )^^|^« M 2 , (23) 

for large enough N. Since the bound is even smaller at earlier times this suggests that the dissipative term is 
negligible, in comparison to the mass term, until time t sp . Equivalently, for weak coupling the damping due to the 
thermal environment has a negligible effect on the quasi-particle mass at these early times. 
Even though ((f> a (x)£(x)) ^ 0, 



< V\ZW^&)) (24) 
(with no summation over a) which, from the behaviour of the noise kernel N, can be bounded as 

(<p a (x)£(x)) £ < fi 2 (j) a (x). 
It follows that </> a satisfies the classical equation 

a<f> a (x) - tf^x) + A0 a (x)^(x) = 0, (25) 

to a good approximation for times t < t sp . For such early times the non- linear term in ( p5| ) can be neglected. For 
times later than t sp neither perturbation theory nor more general non-Gaussian methods are valid. For example, in 
quantum mechanical models (without environment) seventh order calculations in a self-consistent (^-expansion (that 
effectively plays the role of a 1/N expansion in large-N calculations ceases to work once the symmetry has been 
fully broken Jig ), For such times we need to solve for classical vortices, rather as we would solve for the evolution of 
vortex tangles in superfluid i He. 

In terms of the radial and angular fields Eq.(^5|) becomes 

□h + [-/i 2 + Ah 2 - {d^d^e)]\i = 

^(h 2 ^) = 0. (26) 

At early times we have seen that 4>\ and 4>2 are independent, and there are no Goldstone modes. However, by the 
time t « t sp , when h 2 ~ rj 2 and is slowly varying, it follows from the second of the equations that the Goldstone 
modes have appeared. More generally, it is through the coupling of the Goldstone and Higgs (h) modes that classical 
defects appear as solutions to (p5|), in a way that was denied at early times. 

We stress that, in deriving Eq.(|25|), we have not had to restrict ourselves to any particular modes in k-space. 
However, although Eq.(|2^) looks to be valid at all times, it is only a sensible equation once field configurations have 
well-defined probabilities associated to them i.e. after decoherence. Here lies a difficulty in that we have seen that to 
depends on wavelength, with shorter wavelengths taking longer to become classical. Although the long-wavelengths 
with most power in their fluctuations, that determine the separation of vortices, have become classical by time f sp , 
this is an ensemble statement that cannot be applied easily to individual vortex solutions to (|2^). Further, Eq.([2~5|) 
couples <p < to <j) > . The effect of including the short-wavelength modes 0> in the environment is, most likely, only a 
qualitative change at a consistent one- loop level, by augmenting K in ( pl| ) with a comparable term K^. There will be 
a further comparable non-diagonal term \ 2 4>\, J </> a </>b- We anticipate that the Langevin equations for the stochastic 
4> < will be just as one would guess from a mode decomposition of ( |26| ) with short-wavelength modes discarded, and 
terms with K and N equally ignorable. This will be pursued elsewhere. 
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V. CLASSICAL DEFECTS 



A. Line densities and the Gaussian approximation 



As we observed, the most visible signal that the transition has occurred will be the presence of classical topological 
vortices (strings), radial solutions to the equations of motion SS[<p,x] = 0, with \n = 0. These are line defects 
in the field, around which the field phase 6 (<f> — he 10 ) changes by 2ir. We can therefore use line zeroes to track 
classical vortices Plflp"9| . [If we had broken an 0(3) symmetry (a = 1, 2, 3) the corresponding defects would be global 
monopoles, giving qualitatively similar conclusions.] 

The total line-zero density z(x) for the long-wavelength mode fields is |l5| ] 

zj(x) = 5 2 [ ( /» < (x)]|e ufc a^ <1 (x)a fc( /) <2 (x)|. (27) 

With line zeroes as the intersections of planes of zeroes of 4>i and 4>2 separately, we have a phase change of ±2ir 
around the lines of intersection, but for exceptional circumstances. 

^From Eqs. ( [i"4| ) to ( |l6| ) the evolution of the two-field correlator (G< ' (x, 0, t) = G<(r, i)) can be written in terms 
of the unitary (reduced) density matrix p". Thus, 

G<{r,t) = (</><(x)0* (0)) t = ^«6<0<.(x)^< b (O))t 

is well understood. The unstable long- wavelength modes grow exponentially fast. If the power spectrum of the field 
fluctuations is defined by 



d 3 k 

\3 



G<(r,t) = I ----P(k,t)e^ (28) 

Jk<ii { Z7T ) 

then k 2 P(k,t) rapidly develops a 'Bra gg' p eak at k 2 = k\ = (D(n/t). 

While the Gaussian approximation |13j]l4[| is satisfied, the line-zero ensemble density n zelo (t) is determined com- 
pletely by the short- distance behaviour of G< (r, t) as 

«-*>-<««».- ^g|f- (29> 

Since G<(r, t) has short-wavelength modes removed, G<(0,£) is finite. 

Further, for this period when the self-consistent approximation is valid, the field energy (E) t of the system field 4> < 
in a box of volume V becomes 

(i^VKlV^I^ + A^-^l 2 ),) 2 ] 

= V[27m zcro (;)G<(0, t) + \{if - G<(0, t)f] 

= 27rL ze ro(*)G<(0,t) + V[\{if - G<(0,<)) 2 ]. (30) 

As before, x-field fluctuations are absorbed in the definition of /A Eq.(|30|) is obtained by using (|29|), and L zcm (t) = 
Vn zero is the total length of line zeroes, on a scale /i -1 , in the box of volume V. 

We understand Eq. (^p|) as follows. Suppose it were valid from time t = 0, when G<(0, 0) = 0(/i 2 ), until time 
t = t sp . At early times most of the system field energy (proportional to V) is in fluctuations not associated with line 
zeroes, arising from the field potential. As time passes their energy density decreases as the system field approaches 
its post-transition value, becoming approximately zero. In part, this is compensated by the term, arising from the 
field gradients, proportional to the length L zcro of line zeroes, whose energy per unit length increases from 0(/i 2 ) to 
0{rj 2 ). At time t sp , when the fluctuation energy can be ignored, we find |20[] 

(E) t ~ L zoro< 7, (31) 

essentially the energy required to produce a vortex tangle of length L zoro (up to O(l) factors from the logarithmic 
tails). Although these line zeroes have the topological charge and energy of vortices, they are not yet fully-fledged 
defects. Initially, G<(r, t) is very dependent on the value of the cut-off. As a result, line-zeroes are extremely fractal, 
with a separation proportional to the scale at which they are viewed, and are certainly not candidates for defects. 
Once the 'Bragg peak' at k = ko is firmly in the interval k < [i, n zcm (t) becomes insensitive to a cut-off 0(p _1 ). This 
means that line-zeroes are straight at this scale, although they can be approximately random walks at much larger 
scales. For sufficiently weak couplings their density this has happened by time < sp . The final coupling of radial to 
angular modes that turns these proto-vortices into vortices incurs no significant energy change, and n zoro of ( |2^ ) is a 
reliable guide for the initial vortex density. 
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VI. FINAL COMMENTS 



The mechanism for vortex production that we have proposed here has two parts. Firstly, the environment renders 
the order-parameter field classical at early times tr> ^,t sp , by or before the transition is complete. Secondly, classical 
defects evolve from line-zeroes whose resultant density can already be inferred in the linear regime, but whose specific 
attributes are a consequence of the non- linear Langevin equations at the spinodal time t sp . 

This is very different from traditional explanations, and we conclude with a brief summary of them. 

An early explanation, originally due to Kibble |JjJ, and still of common currency, is that thermal fluctuations in 
the Ginzburg regime might lead to the production of vortices, again at early times. The reason is the following: 
once we are below T c , the Ginzburg temperature Tq < T c signals the temperature above which there is a significant 
probability for thermal fluctuations between one degenerate groundstate and another on the scale of the correlation 
length at that temperature. That is, the thermal energy in such a fluctuation matches the energy required to pass 
over the hump of the unstable minimum. This picture presupposes a slow quench, and cannot be accommodated in 
the instantaneous quench approximation that we have used here. However, our suggestion that defects only appear 
at, or about, the spinodal time t sp at a density given by the density of line zeroes is totally at variance with this 
picture, and thermal activation is not the relevant mechanism. 

In fact, this was recognised early by Kibble himself ||, with his later emphasis on strong causal bounds. It had been 
noted earlier Q that the field must behave independently in initially space-like separated regions. When these domains 
with different vacua become causally connected we expect defects to link them. This gives late-time predictions that 
we cannot address with our early-time analysis. A more powerful variant on this theme identifies the time at which 
defects first appear as the time the adiabatic (long-distance) correlation decreases at the speed of light. This time, 
and the resultant density, depends on the quench time tq. There are also difficulties with this in that we have seen 
that the separation of line-zeroes are obtained from the short-distance behaviour of G< (r, t) . The argument was not 
posed for instantaneous quenches of the type discussed above, which correspond to tq = 0(/i _1 ) for which, if taken 
literally, it would give defect formation from time 0(/i _1 ), again incorrect. We are in the process of extending the 
analysis above to slower quenches to make the comparison with the predictions of |S] more useful. 

Subsequently, two approaches have been adopted. In the first, motivated by condensed matter physics, for which 
similar causal arguments apply Q , phenomenological classical stochastic Langevin equations have been proposed (e.g. 
sec [^1|) for the evolution of the fields from early times. These are the counterpart of the phenomenological time- 
dependent Ginzburg-Landau equations of condensed matter and assume classical probabilities from the beginning. 
However, unlike the intermediate-time Langevin equations derived previously, these equations do not have multiplica- 
tive dissipation and multiplicative noise |^ ]l6| , |l7t . The second approach is intrinsically quantum mechanical, treating 
the quantum field as a closed system. Although classical correlations are present in the localisation of the Wigner 
functional p2fl , even for such closed systems, it is difficult to identify defects easily, given that there are no classical 
probabilities at early times. At best, there is late-time dephasing |10|. The mechanism proposed in this paper is 
totally different. 
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